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Abstract 

A simple procedure to estimate ©(af ) and O(a^) corrections to mass-dependent 
observables is conjectured. The method is tested in a number of cases where the 
0{a^) contribution is exactly known, and reasonable agreement is found. 
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Given the great difficulty of evaluating higher-order corrections, it is very desirable to 
have reasonable methods to estimate their sign and magnitude. In fact, significant and 
interesting investigations in this subject have been carried out in the past . 
The aim of this note is to propose a simple estimation method to treat an important class 
of mass-dependent observables. 

We first recall salient features of the estimation methods proposed in the literature. 
Calling R{s) an observable depending on a single time-like kinematic variable s, such as 
a squared centre-of-mass energy, we consider the QCD expansion 



.'2\„nt ,.2\ 
n=0 

where a{ij?) = a^r^\fj?)/n, fi is the renormalization scale, and is the number of 
fiavours active at that scale. By factoring out an appropriate power of s, it is always 
possible to render R{s) dimensionless. Henceforth, we shall adopt this convention. As 
R{s) is renormalization-group invariant, we may choose /i^ = s, in which case Eq. (|l|) 
becomes 

oo 

R{s) = Y.rna''is), (2) 

n=0 

where r„ = Rn{s, s). If R{s) does not depend on masses or other kinematical variables, the 
r„ are numerical constants. Estimations of higher-order corrections using optimization 
procedures based on the fastest apparent convergence (FAC) and the principle of 
minimal sensitivity (PMS) have been carried out in two main scenarios: 
(i) Suppose that tq, ri, and r2 are known. Then the FAC and PMS approaches lead to 
the estimates []T|,3 



PMS 




(3) 



where 

/J„4(ll4,v). A^i(l02-f,v) (4) 

are the first two coefficients of the QCD j3 function. 

(ii) Suppose that ro, ri, r2, and are known. Then in both the FAC and PMS methods 
one finds 



FAC/PMS I n''^3 o'^S '"sA , ^2 \ /-n 

r4 = ra 3 2— - — — + -5- , 5 



where, in the MS scheme, [0] 
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If ra is not known, we may consider employing the estimates from Eq. (|^). Inserting these 
expressions into Eq. one finds 



FAC _ ^ hr^Pi p2 



„PMS 



Formally, the FAC and PMS procedures could be applied in both Minkowskian and 
Euclidean spaces. However, on general grounds it is expected that the optimization 
procedures are more accurate when applied in Euclidean space, as one avoids the presence 
of physical thresholds. Accordingly, when the theoretical expansion for an observable is 
given Minkowskian space, it has been proposed ^ to apply the optimization procedures 
to an associated function defined in the Euclidean region, namely 
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with Q"^ > 0. R{s) admits the inverse integral representation 

1 r-s+ie n(s') 

where D{s') is the analytic continuation of D{Q'^) to the complex plane. Inserting the 
expansion of Eq. into Eq. (H), carrying out the integration, and then setting /x^ = Q^, 
one finds 

oo 

^(Q') = E^^na"(g'), (10) 

n=0 

where 

do = To, di = ri, d2 = r2, (11) 
4 = rs + Y^iPl d4 = u + 7i^po(^r2po + ^npiy (12) 

In the Euclidean approach, one estimates ds and d^ using the expressions analogous to 
Eqs. (H), (I), and (|^) with r„ replaced by dn, and obtains and r4 from Eq. ([T^) . 
At this point, we turn our attention to mass-dependent observables of the form 

oo 

T(.)=m^(/i2)E7^n(3,/^V(/^'), (13) 

n=0 

where m(/i^) is a running quark mass. Again, without loss of generality, we may assume 
that the T„(s,/i^) are dimensionless. Setting jj? = s and defining t„ = T„(s, s), we have 

oo 

T{s) = m\s)Y.tna'\s). (14) 

n=0 
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If T{s) I m'^{s) does not depend on masses or other kinematic variables, the tn are numerical 
constants. The associated function defined in the Euclidean region is 



/•CO 

JO 



Tis) 



s + Q 



2^2 ' 



(15) 



In order to carry out the s integration in Eq. ([TsD, we substitute in Eq. ([T^ ) the expan- 
sions H 



als) 



'PI 



+ 0{a 



(16) 



m{s) = m{fi'') < 1 - a(/i^)7o£ + a\iJ,^)i 



|(/?o + 7o)^-7i 
/5i7o 



+ 
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-y(/5o + 7o) (/3o + I) + + 7i(/5o + 7o)j ^ - 72 

>o + 7o)(/3o + |)(^o + |)^^ 
^/9o + 7o)+7i(/?o + 7o) (/^o + y))^' 

+ 71 (a + y) + 72 (^/?o + 7o)) ^ - 73 + 0(a^^^)| , 



(17) 



where ^ = ln(s//i^) and, in the MS scheme, the coefficients of the mass anomalous dimen- 
sion are [HI 
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Here, ( is Riemann's zeta function, with value ({3) ~ 1.20206. The coefficient 73 is 
presently unknown. For completeness, we have also presented the 0{a^) term in Eq. (|16|), 
although we shall not need it here. These substitutions bring Eq. ([I^) into the form of 
Eq. (P!B| ) with coefficient functions T„(s,/i^) that depend on s only through powers of i. 
Inserting the expression thus obtained into Eq. ([15D and using the elementary integrals 



ds 



{1; 



9. p2. f3. f4, 



} 



TT 



{s + Q 



2\2 
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111 
15 



(19) 



where L = ln((5 //i ), one obtains an expansion of the form 

F(g2) = m2(/i2) ^ F„(Q^/iV(/i')• 

n=0 



(20) 
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Setting ^"^ = Q^, Eq. ( pO] ) becomes 



F{Q')=m'{Q')Y.fna-iQ'), 



(21) 



n=0 

where = Fn{Q'^,Q'^) are numerical constants. Specifically, one finds 

r2 



/o 



TT 



^0, /i=^i, /2 = ^2 + yio7o(/5o + 270), 



h = t3 + y{ti(/5o + 7o)(/5o + 27o)+to[/5i7o + 27i(/5o + 27o)]}, 



(22) 
(23) 



/4 



t4 + vr' |t2(/9o + 7o) [po + ^7o) + h Pi (^/?o + 7o) + ^niPo + 7o) 



+ to 



^270 ^2 ^ 4 

+ 27i(Pi + 7i) + 72 ( Po + ^70 



+ ^to7o(/5o + 7o)(/3o + 270) (y + y 



(24) 



The case of Eqs. (|T|) and (0), in which the m?{^'^) factor is not present, can be obtained 
from Eqs. (p2l)-(p^) by setting 7^ = (2 = 0, 1,2). In fact, the relations between the 
and tn then become identical to those between the dn and r„ in Eqs. (|ll]) and (12). 

Our proposal is to apply the estimation procedure described before to the /„ expansion 
of Eq. ( pTD and to obtain the corresponding t„ coefficients via Eqs. (p3|) and (p^). For 
instance, if the t„ are known for < 2, we obtain /„ for n < 2 from Eq. ( p2|) and estimate 
/s using the expressions analogous to Eq. (|) with r„ replaced by /„. The estimate for 
^3 is then obtained from Eq. (^3]). If the t„ are known for n < 3, we obtain /„ for n < 3 
from Eqs. ( P^D and (P5D, /4 is estimated from the expression analogous to Eq (^ with r„ 
replaced by /„, and t4 follows from Eq. (|2^. If ^3 is not known, we may also attempt to 
estimate from the expressions analogous to Eq. (|^), and once more from Eq. (P^. 
This proposal essentially relates the estimation of the higher-order coefficients in the mass- 
dependent expansion to the previously considered mass-independent case. It should be 
pointed out that there is an element of arbitrariness in this approach. In Eq. (pO]), we 
have set /x' = and proposed to apply the optimization procedure to Z^J^Lo fnoP'iQ'^), the 
cofactor of rn^iQ"^). Had we chosen a different scale /x' 7^ in Eq. (|20|), the expansion in 
Eq. (p^) and the estimation procedure would be different. On the other hand, the choice 
/i' = seems natural and convenient, as all the logarithms in the F„(Q',/i') vanish. 
In fact, this feature has an additional very useful property: it renders the analysis of 
independent of the unknown coefficient 73. 



An interesting application is the estimation of the 0(af) and 0{af) coefficients in 
the evaluation of the partial width T{H — > hadrons) involving final-state quarks with 
running mass m(/i') <^ Mh- The relevant expansion flOl is of the form of Eq. (14) with 
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J- exact 
''3 


j-FAC 
''3 


j-PMS 
''3 


,FAC/PMS 


3 


89.156 


75.729 


80.206 


-945.28 


4 


65.198 


64.956 


68.316 


-1098.8 


5 


41.758 


53.295 


55.547 


-1237.4 



Table 1: Estimations of ^3 and in r{H qq) based on the FAC and PMS optimizations 
of the associated function F{Q^)/m^{Q^), defined in the Euchdean region, and Eqs. (|2^) - 
The estimation of employs the exact value of ^3. 



^/ 


^cxact 
''3 


j-FAC 
''3 


J.PMS 
^-3 


,FAC/PMS 


3 


89.156 


235.82 


240.30 


-527.81 


4 


65.198 


211.20 


214.56 


-748.62 


5 


41.758 


186.67 


188.92 


-949.39 



Table 2: As in Table |T|, but using the original function T{s)/m?{s), defined in the 
Minkowskian region. 

a/s = Mh, to = I, ti = 17/3, and ^2 ~ 35.93996 — 1.35865 n/, where nj is the number 
of active flavours at scale ^/s = Mh- The calculation assumes that there is one massive 
flavour, identical with that present in the final state of the reaction H qq, and nj — 1 
massless ones. Table |I] compares the ^3 estimates of the proposed procedure, based on the 
FAC and PMS optimizations of F(Q^)/m^((5^), with the exact result and also provides 
the ^4 estimate obtained using the exact value of t^. Table ^ displays the corresponding 
estimations based on the optimization of the original function T{s)/rn?'{s) defined in the 
Minkowskian domain. We see that the predicted signs for are correct, but it is apparent 
that the magnitude of the estimations is much closer to the exact result in the Euclidean 
approach. In fact, the estimations of ^3 in Table ^ are fairly good: t\^^ shows errors 
of (—15, —0.4, +28)% for n/ = 3,4,5, respectively; for tg^^ the corresponding errors 
are (—10, +5, +33)%. We also see that the coefficients are predicted to be large and 
negative. 



There are two other cases in which ^3 is exactly known: these are the terms proportional 
to in the absorptive parts of the axial- vector correlators pertinent to the parton-level 
decays — > cs and Z ^ bb |TT|]. Here, rric is evaluated with n/ = 4 at ^/s = Mw and 



with ?7,/ = 5 at ^/s = Mz- In either case, it is assumed that the remaining nf — 1 quarks 
are massless. The corresponding coefficients are ti = 5/3, t2 ~ —3.06004 — 0.02532 ri/ 
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-/.exact 
''3 


''3 


j-PMS 
''3 


,FAC/PMS 


3 


-87.394 


-105.06 


-103.75 


-582.03 


4 


-83.356 


-98.597 


-97.609 


-487.93 


5 


-79.598 


-92.475 


-91.813 


-400.78 



Table 3: Estimations of and in the contribution to r(iy"'' cs) based on the 
FAC and PMS optimizations of the associated function F{Q'^)/m^{Q'^), defined in the 
Euchdean region, and Eqs. (E2D-(|3). The estimation of ^4 employs the exact value of ^3. 



^/ 


j.exact 
^-3 


+FAC 
''3 


+PMS 
^3 


,FAC/PMS 


3 


32.096 


8.6895 


11.587 


-331.80 


4 


20.311 


4.2754 


6.4494 


-400.40 


5 


8.6525 


-0.70475 


0.75256 


-464.41 



Table 4: As in Table ^ but for Z — > bb. 



for W+ cs and h = 11/3, ~ 18.02329 - 0.74754 ri/ for Z bb. Tables | and | 
show the FAC and PMS estimations for these transitions in the Euclidean approach. For 
illustration, we also consider other values of Uf. In the case of —>■ cs, we see once 
more that the estimations are fairly good, with relative errors of 20% or below. Instead, 
in the case of Z ^ bb, the relative errors are large. We note, however, that in this case 
both the exact and estimated ^3 values are relatively small. In fact, a simple, way to 
characterize Tables |I], ^ and ^ is to say that the ^3 estimations have absolute errors of 
order 20 or below. When tf^^'^^ is large, as in Tables |l| and ||, this leads to fairly accurate 
results. 



There are some important cases in which the coefficients are not known. Examples 
include mass relations of the form 

00 

Mg=fi,Y.tna''{fil), (25) 

n=0 

where Mg is the pole mass and fig = mg{fj,g) is the MS mass of quark q. In contrast to 
the previous applications, Eq. ( ^5]) does not depend on an external mass parameter such 
as Mh, M\y, or Mz which, in principle, can have arbitrary values independent of m^. 
On the other hand, defining T(s) = mq{s) Yl'^=Q'tn(^^{s) for arbitrary s > 0, we have the 
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^/ 


j-FAC 
''3 


j-PMS 
''3 


j-FAC 


j-PMS 


3 


152.71 


153.76 


2083.8 


2123.4 


4 


124.10 


124.89 


1544.1 


1571.4 


5 


97.729 


98.259 


1091.0 


1107.8 


6 


73.616 


73.903 


718.74 


727.00 



Table 5: Estimations of and t4 in Eq. (|25D based on the FAC and PMS optimizations 
of the associated function F{Q'^)/m{Q'^), defined in the Euclidean region, and Eqs. (|27|) - 



mathematical identity 



1 1—lJ.q+ie 

M, = — ds' 

Zm J-^iq-ie JO 



ds- 



ns) 



(26) 



[s + s'y 

In analogy with the previous apphcations, we introduce the function F{Q'^) defined by 
Eq. (|15D in the Euchdean region > 0. Because of their hnear dependence on rUg, the 
relations between the and tn are different from those in Eqs. (^^-(^). We now have 

/i = ti, f2 = t2 + — to7o(/5o + 7o), (27) 
6 

^2 



/o 

h 



tn 




7o)(/?o + |)+to 



7o 



7o 



+ t 



+ 



+ 



to 

III 
60 



/32l0 , ll ( n , 71 



//5o , 70 




7i(/5o + 7o) 



:28) 



+ ^(2/5o + 7o) 



^o7o(/?o 



7o)(/?o + | 



(29) 



Assuming once more that nj ~ 1 quarks are massless, we have in the case of Eq. (p5|) 
tQ = l^ti= 4/3, and ta ~ 14.48476 - 1.04137 rij pf. Table | displays the coefficients 
^3 and ^4 in Eq. ( pSj) estimated by the FAC and PMS optimizations of the associated 
function F {Q"^) / m{Q'^) , using Eqs. (H) and (^ with r„ replaced by /„. The entry for 
Uf = 6 corresponding to Mt/ fit, where t is the top quark, gives an expansion very close to 
that obtained in a different approach based on an optimization of the ratios Mt/mt{M^) 



and mt{Ml)liJ,t 



For example, if the expansion is made in powers of a{M^ ) rather 



than a{fif), the FAC result in Table |^ becomes 



1 + ^a{M^) + 8.2366 a^iM^) + 79.838 a^{M^) + 835.69 a^{Mf 



+ 0(a' 



(30) 
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while the approach of Ref. |T3[ leads to 



1 + -a{Ml) + 8.2366 a^(M,2) + 76.172 a^{M'f) + 797.95 a^(M2) 
3 



(31) 



It is well known that, if expressed in terms of /xj, the QCD corrections to Apj, the 



fermionic contributions to the electroweak p parameter, are of the form WA 



Apf 



1 - 0.19325 a{M^) - 3.9696 a^{M^) + O(a^) 



(32) 



where is Fermi's constant. Most of the second-order coefficient in Eq. (p2D, —4.2072, 
arises from the opening of a new channel, namely the double-triangle diagram. It is 



clear that at present there is no sufficient information to optimize Eq. (|32| ). However, if 
one makes the reasonable assumption that the higher-order terms in Eq. (^2]) follow the 
pattern of rather small coefficients shown by the leading contributions, we can combine 
this result with Eqs. ( pOD or ( [3T| ) to estimate and ^4 in Apf/M^. Substituting Eq. (|30|) 
into Eq. (^), one finds 



Apf 



3GjMl_ 
87r2v^ 



1 - 2.8599 a(Mj2) - U.594:a^{M^) - 90.527 a^(M^^ 



- 924.88 a^(M2) + C(a5 



(33) 



The pattern of rapidly increasing coefficients of the same sign displayed in Eqs. (|30D, 
(^), and (^) also emerges from the analysis of infrared- renormalon contributions ||15| . 
As pointed out in Ref. [Q, expansions with much better convergence properties can be 
obtained by optimizing the scale p at which a{p^) is evaluated. 

In summary, we have presented a procedure to estimate C?(a;^) and 0{a\) corrections 
to a class of mass- dependent observables of the type shown in Eqs. ([I3| ) and (|T^) . Although 
there are elements of arbitrariness in its construction, the proposed algorithm, based on 
the optimization of associated expansions in the Euclidean region, is quite simple and 
obviates the dependence on the unknown coefficient 73 of the mass anomalous dimension. 
In the cases where the ^3 coefficients are exactly known, the proposed algorithm estimates 
ts with absolute errors of order 20 or below. In two of the three cases considered, tf^^"^^ 
is large, and the estimations are fairly accurate, with reasonable relative errors. We 
have then generalized the estimation algorithm to important expansions of the form of 
Eq. (^), where the ^3 coefficients are so far unknown. The corresponding ^3 and ^4 
estimations turn out to be close to those found in recent analyses based on alternative 
optimizations procedures. 
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